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ABSTRACT 

The survival of unbound density substructure against orbital mixing imposes strong 
constraints on the slope of the underlying gravitational potential and provides a new 
test on modified gravities. Here we investigate whether the interpretation that the 
stellar clump in Ursa Minor (UMi) dwarf spheroidal galaxy is a 'dynamical fossil' is 
consistent with Modified Newtonian dynamics (MOND). For UMi mass models in- 
ferred by fitting the velocity dispersion profile, the stellar clump around the second 
peak of UMi is erased very rapidly, within 1.25 Gyr (6.5 orbits), even with the in- 
clusion of self-gravity. We find that the clump can hardly survive for more than 2 
Gyr even under more generous conditions. Alternative scenarios which could lead to 
a kinematically cold clump are discussed but, so far, none of them were found to be 
fully satisfactory. Our conclusion is that the cold clump in UMi poses a challenge for 
both ACDM and MOND. 

Key words: galaxies: individual (Ursa Minor dSph) - galaxies: kinematics and dy- 
namics - dark matter - gravitation - stellar dynamics 



1 INTRODUCTION 

The standard concordance cosmological model with cold 
dark matter (ACDM model) is remarkably successful on 
scales larger than 1 Mpc, but it faces challenges on smaller 
scales. For instance, it seems that the theory predicts a too 
cuspy density profile for the dark matter at the centres of 
galaxies (e.g., Trachternach et al. 2008). The Modified New- 
tonian Dynamics (MOND) proposed by Milgrom (1983) has 
proven to be successful in reproducing the kinematics of spi- 
ral galaxies without any assumption of unseen matter (see 
Sanders & McGaugh 2002, for a review), from extremely low 
mass galaxies of low surface brightness (Milgrom & Sanders 
2006) to high luminosity galaxies (Sanders & Noordermeer 
2007) . Gentile et al. (2007) found that the observed rotation 
curves in tidal dwarf galaxies are quite naturally explained 
without any free parameters within MOND, and are incon- 
sistent with the current ACDM theory (see also Kroupa et 
al. 2005) . If MOND is able to naturally account for all of the 
discrepancies faced by ACDM on small scales, this would 
lend strong support to MOND. 

In the Newtonian dark matter scenario, dwarf 
spheroidal galaxies (dSph's) require the largest mass-to-light 
ratios. Hence, dSph's provide a unique testing ground for 
the nature of dark matter and its alternatives (Gerhard & 
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Spergel 1992; Milgrom 1995; Lokas et al. 2006; Sanchez- 
Salcedo et al. 2006; Angus 2008). However, the mass-to-light 
ratios inferred in MOND are very sensitive to uncertainties 
on the structural parameters, luminosities, distances and in- 
ternal velocity dispersions. It is, therefore, important to ex- 
plore other gravitational effects, which, in principle, may 
offer independent tests. 

In the dark matter paradigm, the velocity dispersion 
profiles of the brightest dSphs are consistent with both a 
cuspy NFW halo and a cored dark halo. However, there ex- 
ists some indirect evidence that dSph galaxies may possess a 
core (Kleyna et al. 2003, hereafter K03; Goerdt et al. 2006; 
Sanchez-Salcedo et al. 2006). In particular, K03 considered 
the survival of the cold density substructure detected in pho- 
tometric data in Ursa Minor (UMi). This substructure ap- 
pears as an off-centre localized stellar clump with low ve- 
locity dispersion ~ 0.5 km s _1 . K03 concluded that the sec- 
ondary peak in UMi is a long-lived structure, surviving in 
phase-space because the underlying gravitational potential 
is close to harmonic. In the standard dark matter paradigm, 
this implies that the dark matter halo in UMi must have 
a large core because, if the dark halo has a central density 
cusp, the clump should have diluted in ~ 1 Gyr. Even if 
dSph galaxies are influenced to some degree by the tidal 
forces exerted by the Milky Way, it is unlikely that such a 
large core may have a tidal origin (e.g., Stoehr et al. 2002; 
Hayashi et al. 2003; Penarrubia et al. 2008). It is worthwhile 



Sanchez- Salcedo & Lora 



exploring if the competing MOND scenario can explain the 
survival of cold substructures in dSph galaxies in a natural 
way. The question that arises is whether the MOND gravita- 
tional potential can mimic the potential of a dark halo with 
a core in order to explain the very longevity of the dynami- 
cal fossil in UMi. In other words, can the inference of cored 
haloes in dSphs be accomodated naturally within MOND? 
This paper is organized as follows. In section 2, we give 
a general statement of the problem and outline some key 
analytic results for tidal dissolution of satellite systems. In 
section 3, we briefly describe the observational properties of 
UMi dwarf and its clump. Some important issues regard- 
ing the assumptions and approximations made to study the 
evolution of the clump in MOND framework are given in 
section 4. Section 5 presents results on the evolution of a 
stellar clump in MOND. In section 6, we discuss alternative 
scenarios to account for the persistence of kinematically cold 
substructure in dSph galaxies. Finally, our conclusions are 
summarized in section 7. 



2 EVOLUTION OF SATELLITE SYSTEMS 

AND UNBOUND CLUMPS: TIDAL RADIUS 
AND EPICYCLIC THEORY 

Suppose that a small clump of mass M c i is on a circular orbit 
with radius 7? c i r around its host galaxy, which is assumed to 
be spherical. If the circular speed with which a test particle 
orbits in the gravitational potential of the host galaxy can 
be described by a power-law v% oc r l with — 1 < I < 2, the 
tidal radius is given by 



r t — k 



Xc 



-,1/3 



Rcir 



(1) 



where Xcir is the ratio between the mass of the clump and the 
mass of the host galaxy mass inside Rdr, and k is Keenan's 
(1981a, b) factor, which accounts for the elongation of the 
zero- velocity surface along the line joining the clump and the 
centre of the host galaxy. For the Keenan factor, we will take 
the nominal value k — 1 (see section 2 in Sanchez- Salcedo & 
Hernandez 2007 for a discussion). Zhao (2005) and Zhao & 
Tian (2006) showed that the expression for the dynamical 
tidal radius given in Eq. (1) is also valid in MOND as long as 
Xcir is taken as the ratio between the "true" baryonic masses 
of the satellite and the host galaxy (see Sanchez-Salcedo & 
Hernandez 2007 for a discussion). We see from Eq. (1) that it 
may happen that the clump is essentially unbound (M c i — > 
and Xcir — > 0), but it is not dissolved because I —¥l. This is 
true in both Newtonian and MOND gravity. The case I = 2 
corresponds to the harmonic potential. If the potential is 
quasi-harmonic, the orbital period depends weakly on the 
energy of the stars and hence the phase mixing is highly 
suppressed even if self-gravity of the clump is ignored. 

Much of early work on the dynamics of unbound satel- 
lites used the epicycle theory to describe their evolution 
(e.g., Oort 1965; Wielen 1977; Fellhauer & Heggie 2005). 
From the linear epicyclic theory, we learn that it is possible 
to prepare a small unbound system with a nonzero velocity 
dispersion that will orbit with almost no secular spreading 
(e.g., Kuhn 1993). In fact, in an idealized situation where 
all the stars of the clump have the same guiding centre, 
the clump would preserve its original size in the epicyclic 



approximation because the epicycle and orbital frequencies 
are the same for all stars. However, even if we start with a 
very idealized situation in which all the stars of the clump 
move on retrograde orbits resembling epicycles with a com- 
mon guiding centre, self-gravity of the clump destroys the 
coherent stellar motions and the dissolution time is infinite 
only when the density of the clump is negligible. Fellhauer 
& Heggie (2005) studied the time of dissolution of such un- 
bound stellar system, in Newtonian gravity, as a function 
of its initial internal density and size. They found that a 
system orbiting within a galaxy modeled as an isothermal 
sphere, can survive for about 10 galactic orbits only if it 
has a very low density (~ 5% of the background density) 
and small radius (0.5% of the orbital radius). For larger in- 
ternal densities of the clump, the dissolution becomes more 
rapid. In the Appendix A, we extend the theoretical anal- 
ysis of Fellhauer & Heggie (2005) to MOND gravity. We 
consider a small unbound clump orbiting on a circular orbit 
around the centre of a galaxy. It is shown that, because of 
the enhanced self-gravity in MOND, the dissolution time of 
a certain unbound clump is shorter in MOND than in its 
equivalent Newtonian galaxy, that is, the Newtonian galaxy 
(the same distribution of stars plus additional dark matter) 
which has the same "dynamics" as the MOND galaxy. 



3 URSA MINOR AND ITS DYNAMICAL 
FOSSIL 

3.1 UMi: Structural parameters 

Ursa Minor is located at a heliocentric distance of D — 76 ± 

4 kpc (Carrera et al. 2002; Bellazzini et al. 2002). The star 
formation history and the characteristics of the population 
of blue straggler stars suggest that UMi is a truly 'fossil' 
galaxy, where star formation ceased completely more than 8 
Gyr ago (Carrera et al. 2002; Mapelli et al. 2007). The shape 
of the inner isodensity contours of the surface density of stars 
appears to be elliptical with a large ellipticity (e = 0.54). 
The King core radius of the stellar component along the 
semimajor axis is 17.9' (~ 395 pc at 76 kpc) (Palma et 
al. 2003). 

In the Newtonian dark matter scenario, values for 
M/Lv ranging between 30 to 200 can be inferred depend- 
ing on the adopted global luminosity, distance and veloc- 
ity dispersion profile. In particular, if we rescale the results 
of Gilmore et al. (2007) for D — 76 kpc and use the up- 
dated luminosity of L v = 1.1 x 10 6 Lq (Palma et al. 2003), 
the total mass within 0.6 kpc is > 7 x 1O 7 M0, resulting in 
a M/Lv > 64Mq/Lq. If we use the dispersion profile of 
Munoz et al. (2005), the M/Lv is reduced by a factor ~ 2. 



3.2 The secondary peak of UMi revisited 

Several authors have reported structure in the surface den- 
sity of UMi on a variety of scales (Olszewski & Aaron- 
son 1985; Demers et al. 1995; Irwin & Hatzidimitriou 1995; 
Kleyna et al. 1998; Eskridge & Schweitzer 2001; Palma et 
al. 2003). The most remarkable feature is the double off- 
centered peak in the stellar surface number density. In this 
paper, we concentrate on the secondary peak, observed on 
the northeastern side of the major axis of UMi. 
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Figure 1. Snapshots of an unbound clump at t = 0, 2, 4.8 and 9 Gyr in a dark matter halo with a scale radius R s equal to 2.4 times 
the semimajor axis of the size of clump's orbit. The clump has an initial la radius of 12 pc, one-dimensional velocity dispersion of 0.5 
km s _1 and a speed at apogalacticon equal to 0.07 times the circular speed. 



K03 found that the radius of the clump is > 6' (equiv- 
alent to 135 pc at D — 76 kpc); this is the largest circle 
where the isophlet contours indicate the presence of a sec- 
ondary peak in the stars. The velocity of stars that form the 
secondary peak is best fitted by a two-Gaussian population, 
one representing the 8.8 km s _1 Gaussian background and 
the other representing a cold subpopulation of 0.5 km s _1 
velocity dispersion. The cold group of stars has the same 
mean velocity as the systematic velocity of UMi, so it is 
likely that its orbit is close to the plane of the sky. 

The luminosity bump can be fitted approximately by a 
Gaussian density profile p — po exp(— r' 2 /2ro), where r' is 
the radial distance to the peak in the stars, and ro — 35 pc 
(e.g., Palma et al. 2003). The projected surface density of 
the clump along its centre is: -\/2irporo- Since the underlying 
stellar population of UMi also contributes to the observed 
surface density, we need the fraction / of stars that belongs 
to the kinematically cold subpopulation. K03 found that the 
subpopulation fraction is / = 0.7 in the best fit. Equating 
V^poro with / times the observed surface density at the 
secondary peak, we can infer po and then the total mass in 
the clump M c i ~ 16po^o- Using / = 0.7, we obtain 



M cl = 7.8 x 1O 4 M ( ^- 



(2) 



where Y* is the V band mass-to-light ratio of UMi stellar 
population. For a normal stellar population Y* « 2, we have 
M cl = 2.5 x 1O 4 M . 

An inspection of the stellar isophlets in UMi suggests 
that the secondary peak is unbound material because of its 
large angular size in the plane of the sky, as compared to 
the expected value for a bound cluster, and its bend in the 
isodensity contours, indicative of tidal disruption (Palma et 
al. 2003). Let us estimate the tidal radius of a stellar clus- 
ter of mass M c i — 2.5 x 10 4 Mq, as UMi's clump, on a 
circular orbit in a galaxy with a flat rotation curve at Rdi 
(i.e. I = in Eq. 1). For a UMi- like galaxy, the mass interior 
to .Rcir = 390 pc in the standard dark matter scenario is 
~ 10 7 M©. According to Eq. (1), the clump would remain 
gravitationally bound only if all the mass that comprises the 
stellar cluster were placed within a sphere of 42 pc radius. 
The appearance of such a hypothetical cluster should re- 
semble a globular cluster as those detected in Fornax dSph 



galaxy and, therefore, easily detectable as a distinct spheri- 
cal bound object. This is not the case for the secondary peak 
in UMi and, thus, it turns out that the clump is unbound 
at least in the dark matter scenario. 

Since there is no evidence that the regions of density ex- 
cess have stellar populations that differ from the main body 
of UMi (Kleyna et al. 1998), the most plausible interpreta- 
tion is that this density structure, or clump, is a disrupted 
stellar cluster (K03; Read et al. 2006; section 6). K03 argues 
that the persistence of a dynamically separate unbound en- 
tity for 10 — 12 Gyr, is possible only if the background po- 
tential is quasi-harmonic (see section 2), implying that the 
dark halo has a large core. For illustration, Figure 1 shows 
the evolution of a clump made of test particles, orbiting in 
a spherical dark halo with a density mass profile of the form 



pdm(r) 



pmax 



(l + [r/R s ] 2 y/2 : 



(3) 



where R a is the scale radius. Following K03, the central den- 
sity p max was selected such that the circular velocity is 19 
km s _1 at an angular distance of 31.5', i.e. at ~ 700 pc. 
The contribution of the stellar mass to the gravitational 
potential was ignored because its mass interior to 0.4 kpc 
(0.9 x 1O 6 M0) is only 9%, or less, of the dark matter mass, 
if we assume that UMi has a normal stellar U-band mass- 
to- light ratio of ~ 2Mq/Lv,q- The clump was dropped at 
a distance R s /2.4 on a near radial orbit with a tangential 
velocity equal to 0.07 times the circular velocity at its initial 
position. The clump has initially a la radius of 12 pc and 
one-dimensional velocity dispersion of 0.5 km s" 1 , similar to 
those selected by K03 in order to facilitate comparison. The 
orbit lies in the (x,y) -plane. From Figure 1, we see that, 
although the density structure becomes more extended over 
time, it may survive for many Gyr if we assume that the 
scale radius of the dark halo is at least 2-3 times the size of 
clump's orbit, confirming K03 results. Due to the phase mix- 
ing, those components of the velocity dispersion along the 
orbital plane (i.e. a x and a y ) increase with time, whereas 
the perpendicular component to the orbital plane remains 
constant. As a consequence, the size of the clump does not 
increase along the axis perpendicular to the orbital plane, 
and remains dynamically cold in that direction (a z ~ 0.5 
kms^ 1 ). 
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Lora et al. (2009) have simulated the dynamics of a 
clump including self-gravity. They show that density sub- 
structure can persist for ~ 12 Gyr in a halo with a scale 
radius equal to 1.5 times the size of the clump's orbit. As- 
suming that the size of the orbit is ~ 390 pc, a scale radius 
for the dark halo of R 3 = 1.5 x 390 pc = 580 pc is required. 
Note that R 3 was defined as the radius at which the density 
has dropped a factor \/2. This large core is at odds with the 
predictions of ACDM. Therefore, the interpretation that the 
clump survives because the gravitational potential is quasi- 
harmonic poses a challenge to the standard ACDM, which 
predicts cuspy NFW profiles. In the next sections, we inves- 
tigate whether MOND can account for the survival of the 
clump. 
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THE CASE OF MOND: 
APPROXIMATIONS 



CONCEPTS AND 



In MOND, the luminous density profile determines the 
shape of the gravitational potential and, thus, the index I in 
Eq. (1). Moreover, Xcir is not sensitive to the adopted stel- 
lar mass-to-light ratio as soon as we assume that the stellar 
mass-to-light ratio of the clump is similar to that of the stel- 
lar bulk of UMi. Consequently and according to Eq. (1), the 
tidal radius is very robust to the adopted stellar mass-to- 
light ratio. Since the only free parameter is, in principle, the 
stellar mass-to-light ratio, the survival of the clump in UMi 
is an interesting dynamical test to MOND. We must warn 
here that Eq. (1) was derived for MOND satellite systems 
embedded in isolated host galaxies. As it will become clear 
later, UMi is affected by the Galactic field and cannot be 
treated as it was in isolation. 

In MOND framework, the gravitational potential de- 
scribing the force acting on a star in UMi obeys the modified 
Poisson equation of Bekenstein & Milgrom (1984) 



V ■ {p(x)V$) = AwGp, 



(4) 



where x — |V$|/ao, do — 1-2 x 10~ 8 cm s~ 2 is the uni- 
versal acceleration constant of the theory and p(x) is the 
interpolating function which runs smoothly from j-i(x) = x 
at x <C 1 to p(x) = 1 at x 3> 1. Equation (4) must be 
solved with the boundary condition V$ — > — g E , where g E 
is the external gravity acting on UMi, which has a magnitude 
gE = V 2 /R gc , where V is the Galactic rotational velocity 
at R gc , which coincides with the asymptotic rotation veloc- 
ity Uoo for the Milky Way. Although Uoo is very difficult to 
determine observationally in our Galaxy, Famaey & Binney 
(2005) and Sanchez- Salcedo & Hernandez (2007) argue that 
it must be of 170 ± 5 km s~ J by adopting a mass model for 
the Milky Way under MOND. 

A clump star feels the external acceleration created by 
the Milky Way (g E ), the acceleration generated by UMi 
background stars (g z ) and the acceleration generated by all 
the other stars of the clump (g int ). Due to the non-linearity 
of the MOND field equation, the gravitational acceleration 
g t is altered by the gravitational field of the Milky Way, 
whereas g int is altered by both g E and g x . 



Figure 2. The observed UMi line-of-sight velocity dispersion 
versus projected radius is shown as diamonds with errorbars. 
The solid line corresponds to the best fitting model in MOND 
(Y* = 5.8) and the dashed line shows the predicted velocity dis- 
persion for Y* = 2.2, which is at la from the best fitting value. 



4.1 Modelling UMi background gravitational 
potential under MOND 

Only for very special configurations (one-dimensional sym- 
metry -spherical, cylindrical or plane symmetric systems- 
or Kuzmin discs) the MOND acceleration g is related to 
the Newtonian acceleration, g N , by the algebraic relation 
p(\g\/a )g = g N (Brada & Milgrom 1995). Milgrom (1986) 
showed that when the mass distribution involves a small 
perturbation in a medium under a uniform external field, 
the perturbation in the gravitational potential satisfies an 
anisotropic Poisson-like equation with a dilation along z, 
which is taken as the direction of the external field g E 
(see the Appendix A). For instance, the 'modified' Plum- 
mer model for a satellite of mass M in this quasi-Newtonian 
limit is (e.g., Sanchez-Salcedo 2009): 



p P (r) = p c 



x + y 

1 + 5^- + 

rp 



(1 + Lo)r| 

GM 



£*o \/(l + L )(x 2 + y 2 ) + z 2 + f 2 



where 



Pc 



3M 



4tt(1 + L y/ 2 rl : 



(1 + L 



,V2, 



-5/2 

(5) 
(6) 

(7) 
are the characteristic radii and 



Lq is the logarithmic derivative of po = p(gE/ao) and hence 
< Lq < 1. This pair is solution of the hydrostatic Jeans 
equation for an isotropic velocity dispersion tensor. We see 
that the potential becomes Newtonian, but with a larger 
effective gravitational constant p^G, and anisotropic. In 
the case of halo objects at distances D much larger than the 
galactocentric distance of the Sun ~ 8 kpc, the direction of 
the external field from the Milky Way almost coincides with 
the line-of-sight. 

As an approximation and effective way to take into ac- 
count the external field effect (EFE) in dSph galaxies, Angus 
(2008) used the following relation between the internal g T 
and external g E accelerations: 
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SiMz) = 9am 



(8) 



where g N is the Newtonian internal acceleration of the sub- 
system and x = \J g\ + g%/ao- In view of the success of the 
model in reproducing the kinematics of the classical dSph 
galaxies (Angus 2008), we will use Eq. (8) to derive g 7 . As we 
will see in the following, Equation (8) is interpreted as one of 
the simplest way to interpolate between the isolated and the 
quasi-Newtonian regimes. Due to its simplicity, it has proven 
to be very useful to treat the EFE (e.g., Wu et al. 2007; 
Angus & McGaugh 2008; see also Llinares et al. 2008). At 
locations where gi 3> gE, X — gi/ao and Eq. (8) can be 
simplified to the algebraic relation n(gi/ao)g I = g N , which 
is the exact solution of the Bekenstein-Milgrom equation if 
the subsystem is spherical or cylindrical (and it is a good 
approximation in other geometries like axisymmetric expo- 
nential discs, e.g., Brada & Milgrom 1995). In the opposite 
regime, i.e. in the quasi-Newtonian limit, where the dynam- 
ics is dominated by the external field, that is gE 3> gi, 
Eq. (8) can be written as ^{gE/a-ojgj = g N . Thus this 
approximation recovers the quasi-Newtonian dynamics pre- 
dicted in linear theory, but does not capture the anisotropic 
dilation along the external field direction. 

Suppose that we have a stellar clump embedded in an 
ellipsoidal potential, such as those described by Eq. (6). This 
flattened potential creates a torque on our clump, when 
the orbital plane of the clump is misaligned with the in- 
stantaneous direction of the external field (Wu et al. 2007). 
This causes additional orbital mixing that contribute to the 
dissolution of unbound clumps. Since the approximation of 
Eq. (8) cannot capture dilation, it neglects this source of 
orbital mixing. A more refined model to include the EFE 
dilation is not warranted, due to our ignorance about the 
density mass profile of UMi along the line of sight. 

Angus (2008) used Jeans analysis to model the UMi ve- 
locity dispersion profile versus circular radii (rather than el- 
liptical radii) reported in Munoz et al. (2005) under MOND 
and demonstrated that the best-fit model for UMi's veloc- 
ity dispersion is very reasonable. Therefore, we will adopt 
exactly the same parameters for UMi as in Angus (2008), 
which are compiled in our Table 1. The mass density of UMi 
was modelled by a spherical King model with a core ra- 
dius of 17.9', which corresponds to the core radius along the 
semimajor axis. The velocity anisotropy was taken variable 
but the fit turned to be quite compatible with a constant 
anisotropy of 0.7. Figure 2 shows the observed velocity dis- 
persion from Munoz et al. (2005), together with the best 
MOND fit (T* = 5.8, where again T* denotes the stellar 
U-band mass-to-light ratio in solar units), as well as the 
predicted MOND velocity dispersion profile when the la er- 
ror in T* is considered (T* = 2.2) 1 . 

In the best fitting model, the circular velocity of a test 
particle at the stellar core radius (r c = 395 pc) of UMi is ~ 
13 km s _1 (see Fig. 3). Therefore, the characteristic internal 
acceleration, [v c (r c )] 2 /r c ~ 0.14 x 10~ 8 cm s~ 2 , is much 
smaller than ao — 1-2 x 10~ 8 cm s~ 2 . On the other hand, 
the external acceleration, V^/R gc — 0.12 x 10~ 8 cm s~ 2 , is 



The discrepant values for Y* between Angus (2008) and 
Sanchcz-Salccdo & Hernandez (2007) is traced to the different 
distance, luminosity and velocity dispersion profile used (see sec- 
tion 3.1). 




Figure 3. Circular speed curves for different models of UMi. For 
the reference MOND model, which has Y* = 5.8, the solid line 
shows the MOND rotation curve, whereas the short dashed line 
gives the Newtonian circular speed multiplied by a factor 2.5 to 
show that both curves almost match. Also shown is the rotation 
curve for a model with the same parameters than the reference 
model but with Y* = 2.2 (dot-dashed line). The dotted line gives 
the circular speed in the reference case (Y* = 5.8) but when the 
external field strength is doubled. For comparison, the rotation 
curve in Newtonian gravity with the dark matter halo described 
in section 3.2 is also plotted (long dashed line). 



also much smaller than ao- We conclude that UMi internal 
dynamics lies in the deep-MOND regime and thus it is not 
sensitive to the exact form of the interpolating function. The 
ratio between the internal acceleration at UMi's core radius 
r c and the external acceleration is a good measure of the 
importance of EFE. For our reference parameters given in 
Table 1, this ratio is ~ 1.1, which implies that UMi dSph 
galaxy is at an intermediate regime; it can be considered 
neither as isolated nor as external field dominated. 



4.2 Self-gravity of the clump 

We may define the internal acceleration g int of a clump's star 
as g int — — V$ — g E — g x . In fact, the clump is embedded 
in the external field generated by all the other particles of 
UMi and the Milky Way. Due to the EFE, g int depends on 
both accelerations g E and g T . Suppose for a moment that 
the clump was spherical and seated at the centre of UMi. 
Since gi is small compared to gE in the centre of UMi, we 
should recover an equation identical to Eq. (8) for g int , i.e. 



SintMz) = 9in 



(9) 



where x — \/gf nt + g%/ a o and g int N is the Newtonian in- 
ternal acceleration of a star in the clump. If the clump is 
displaced at those distances from the centre of UMi where 
gi is not longer negligible relative to gE, Equation (9) over- 
estimates clump's self-gravity, making it more resilient to 
tidal destruction. 

At g^t S> gE, Equation (9) is the solution of the mod- 
ified Poisson equation when our clump is spherical. There- 
fore, Equation (9) is valid provided that the mass density 
distribution of the clump is close to spherical at g^t 3> gs- 
As time goes by, the clump and its tidal debris will form a 
non-spherical mass configuration because the stars lost from 
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Table 1. A summary of relevant parameters for the reference model. 
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supported by observations —see text—). 



0.4 

0.2 

0.0 

-0.2 

-0.4 

0.4 

0.2 

0.0 

-0.2 

-0.4 



t=o :: 



0.2 Gyr:: 



t I I I I I I I t I I I t 



\ 



1 Gyr :: 




hii l lulu 'I 




\ 



2 Gyr 



• . . :.■■: t *ife ■', ■ 
■.: .■ ■■_■ v.'. ,.' . 

' '■y-lii'^'.'i:-'-. 
*:. , '.J-v--:-U:-i-. < -*, 



t i i i t i i i l i i i l 



/ 



-0.4 -0.2 0.0 0.2 0.4 -0.4 -0.2 0.0 0.2 0.4 -0.4 -0.Z 0.0 0.Z 0.4 -0.4 -0.2 0.0 0.2 0.4 
v [kpc] x [kpc] x [kpc] x [kpc] 



Figure 4. Evolution of a clump in the reference MOND model for a near radial orbit (top panels) and on a circular orbit (bottom panels) 
when clump's self-gravity is ignored. The orbital plane of the clump lies in the (x, y) plane. Times since the start of the simulation are 
indicated. 



the clump will form two tidal tails, one leading the clump 
and one trailing it. However, the break of spherical symme- 
try does not represent a real limitation because, as we will 
see in the next section, gi nt is comparable to qe only at the 
very beginning of the simulations when the clump is roughly 
spherical. In fact, internal accelerations of the order of qe 
only occur at initial times and inside the clump, where the 
mass distribution can be considered as spherical. 

In order to quantify the role of self-gravity, we will 
carry simulations in which the self-gravity is turned off (non- 
selfgravitating case, NSG) and simulations in which self- 
gravity is included through Eq. (9) . According to the ongo- 
ing discussion, these are the limits of what can be expected 



in MOND; the survival times derived with Eq. (9) represent 
upper values. In the remainder of the paper, we will refer to 
Eq. (9) as the reduced EFE approximation (REFE). Direct 
N-body simulations of 1200 particles were carried out using 
the code described in Lora et al. (2009) with the REFE ap- 
proximation. The convergence of the results was tested by 
comparing runs with different softening radii and number of 
particles. 
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Figure 5. Same as Figure 4 but for a MOND model with T* = 2.2. 



5 EVOLUTION OF A STELLAR CLUMP 

5.1 NSG simulations 

We consider first the reference model in which T* = 5.8. 
Figure 4 shows the disintegration of a NSG clump on a near 
radial orbit with an apogalactic distance of 390 pc and on 
a circular orbit with a radius of 390 pc. In the case of a 
near radial orbit, the clump was dropped at apogalacticon 
with a tangential velocity equal to 0.07 the local circular 
velocity. The initial la radius of 12 pc and velocity disper- 
sion of the stars in the clump of 0.5 km s _1 are the same as 
those used in section 3.2. Phase mixing dissolves the clump 
very rapidly; the clump doubles its size in 0.2 Gyr (essen- 
tially one orbital revolution of the clump about UMi) and 
it becomes completely diluted in less than 1.25 Gyr (~ 6.5 
orbital periods). When the clump is on a circular orbit, the 
tidal debris populates the whole ring of the orbit, due to 
the effect of differential rotation. The clump is shortlived 
because the MOND circular velocity v c (r) of UMi is essen- 
tially a scaled-version of the Newtonian (without dark mat- 
ter) circular velocity at r < 0.6 kpc (see Fig. 3). Thus, the 
evolution of the clump in MOND resembles that in a purely 
Newtonian galaxy with a dark halo having a core equal to 
r c , i.e. a mass-follows-light model. Since the core radius of 
the fictitious dark halo is almost equal to the clump's orbit, 
the group of stars is erased in a few crossing times. In other 
words, MOND is not able to change significantly the shape 
of the gravitational potential within the orbit of the clump, 
at least not enough to generate a near-isochrone potential 
needed to guarantee the survival of the clump. 



Since the destruction time roughly scales as the galactic 
crossing time for the clump, one can increase the longevity of 
the clump by decreasing the mass of UMi if a lower mass-to- 
light ratio of the stellar population of UMi is adopted. Figure 
5 shows the evolution of the same clump if the body of the 
stars in UMi has T* = 2.2. The clump is diluted in < 2 Gyr. 
Even when using a M/Ly value seven times smaller than the 
preferred mass-to-light ratio (i.e. T* = 0.83) and the clump 
is initialized with a velocity dispersion of 0.25 km s _1 , the 
group of stars dissolves in < 2.5 Gyr (see Fig. 6). Values 
for the mass-to-light ratio below 0.83 may be problematic 
because then UMi may be too fragile to the process of tidal 
disruption. 

Similarly, the survival time can also be enhanced by as- 
suming a larger value for the external field acceleration §e- 
When we double the strength of the external field in our 
reference model, the clump still dissipates in < 2 Gyr. For 
comparison, the circular speed for this model was also plot- 
ted in Fig. 3. The uncertainty on V x cannot change signifi- 
cantly the disruption timescale. Doubling §e corresponds to 
an asymptotic circular velocity of 250 km s^ 1 , which sounds 
unrealistic in MOND. 



5.2 Simulations with self-gravity 

In the previous section, we have considered the fate of an 
unbound clump in MOND. Even if the clump is initially very 
compact (a la radius of 12 pc), orbital phase mixing dilutes 
the overdensity structure within less than 2 Gyr in MOND 
models. When self-gravity of the clump is included, a clump 
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Figure 6. Distribution of the stellar debris projected onto the 
orbital plane at t = 2.5 Gyr in a MOND case with Y* = 0.83, 
when the clump with initial velocity dispersion of 0.25 km s _1 
is on a near radial orbit (upper panel) and on a circular orbit 
(bottom panel). Self-gravity of the clump was turned off. 



with an initial la radius of 12 pc remains too compact to 
explain its present appearance (see Fig. 7). Therefore, we 
have examined the disruption timescale of a clump with la 
radius of 35 pc, which is similar to the present extent of 
the clump. A compilation of the parameters for our refer- 
ence model is given in Table 1. When clump's self-gravity 
is taken into account, a velocity dispersion of 0.5 km s _1 is 
not enough to keep the clump in virial equilibrium and the 
clump collapses to a hotter compact configuration. In order 
to have the clump in a relaxed state, we needed to adjust 
the initial one-dimensional velocity dispersion. 

Figures 8 and 9 show the evolution of the clump (assum- 
ing T* = 5.8) on a near radial orbit and on circular orbit as 
those described in section 5.1 but including self-gravity in 
the REFE approximation. The initial one-dimensional ve- 
locity dispersion of the clump is 1 km s _1 . Although the 
disruption of the clump is slightly faster when the clump is 
on a near-radial orbit, dissolution of the clump occurs within 
1.25 Gyr in both cases. When self-gravity is turned off, the 
cluster of stars is erased within 0.5 Gyr. 

In order to gain more physical insight, we will focus on 
the simulation with the clump on circular orbit. The clump 
cannot be considered in isolation at any time because the 
mean internal acceleration due to mutual forces between the 
stars of the clump, (<?mt), is initially ~ 0.9gE and decreases 
in time as the stars of clump becomes more loosely bound. 
For instance, (gint) averaged over all the simulated stars is 
0.37g E at t = 0.5 Gyr. In Fig. 10, the MOND boost fac- 
tor gint/flint.N at t — 0.5 Gyr is plotted as a function of the 




Figure 7. Stellar surface density in Mq pc on the orbital 
plane after 6 Gyr, for a clump with initial 1<t radius of 12 pc, 
one-dimensional velocity dispersion of 1.8 km s and a speed at 
apogalacticon equal to 0.07 times the circular speed. Self-gravity 
of the clump was included. 



distance to clump's centre. Once a star has been stripped 
away from the clump, at distances > 100 pc, it feels a con- 
stant boost factor and its dynamics is Newtonian but with 
a larger effective Newton constant G. All the particles lie in 
the quasi-Newtonian regime when the clump is erased. 

Newtonian conservation laws such as the the conserva- 
tion of the total energy and total angular momentum are 
only preserved in the REFE approximation when the dy- 
namics of all the stars is quasi-Newtonian. In order to quan- 
tify the degree of conservation in our simulations, we have 
computed the net internal force m ~^2 N g int (being m the par- 
ticle mass) summed over all the simulated particles and the 
total angular momentum. The net internal force fluctuates 
in time but it is always less than 0.5 percent the character- 
istic background force M c iv^ c /r c , where v cc is the circular 
speed around UMi at r c . The total angular momentum, on 
the other hand, increases by 4 percent during the first Gyr 
and by 0.5 percent between 1 Gyr and 1.5 Gyr. Therefore, 
the inspiraling of the clump from an initial radius of 390 
pc to 355 pc at 1 Gyr (see Fig. 9) is not an artifact but 
the result of the gravitational transfer of angular momen- 
tum from the clump to the trailing tidal tail. We must also 
notice at this point that the density substructures displayed 
in the map at 1.5 Gyr do not have statistical significance, 
as Fig. 11 demonstrates. 

A reduction of Y* leads to an increase of the crossing 
time of the clump around UMi's centre but also the mass 
of the clump is lowered if we assume that the clump has 
the same Y* as the background stars. For Y* = 2.2 (M c / = 
3 x 10 4 M©), the clump dissolves within 1.25 Gyr as well 
because clump's self-gravity is less relevant. 

It may be worthwhile exploring whether the destruction 
timescale depends on the adopted value of the apocenter of 
the clump's orbit. For an apogalactic distance of 750 pc, 
we find that the clump is erased in < 1 Gyr (see Fig. 12). 
Therefore, the persistence of the clump cannot be accounted 
for by altering its orbital parameters. 
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Figure 8. Evolution of the stellar surface density in Mq pc -2 for a clump with initial la radius of 35 pc, one-dimensional velocity 
dispersion of 1 km s — 1 and a speed at apogalacticon equal to 0.07 times the circular speed. Self-gravity of the clump was included. 



Self-consistent simulations of the evolution of a low- 
mass satellite were carried out by Kroupa (1997) and 
Klessen & Kroupa (1998). They found that the unbound 
distinguishable remnant, consisting of particles with phase- 
space characteristics that reduce spreading along the orbit 
as those described by Kuhn (1993), contains a mass of about 
1% of the initial mass of the progenitor. In the case of UMi's 
clump, which has a mass of ~ 7.8 x 1O 4 M0 for T* = 5.8, 
the long-lived remnant should have a mass of ~ 780Mq , too 
small to be detected in our simulations. 

In all our previous REFE models, the external acceler- 
ation was taken as constant and equal to the instantaneous 
external acceleration of UMi at its current position. The 
external acceleration, however, may increase with time as 
UMi orbit penetrates deeper into the Milky Way halo when 
it moves towards pericentre. Piatek et al. (2005) measured 
the proper motion for UMi and found a perigalacticon and 
apogalacticon of 40±3o kpc and 89±J} kpc (at 95% con- 
fidence level), respectively. According to Eq. (8), as UMi 
plunges in on an eccentric orbit, it puffs up and expands adi- 
abatically because gi gets smaller at larger external accelera- 
tions (Brada & Milgrom 2000). Thus, the rotation timescale 
of the clump around UMi slows. For the same reason, the 
self-gravity of the clump lowers and the clump also inflates, 
rendering itself more susceptible to tidal dissolution. In other 
words, when the strength of the external acceleration in- 
creases, the internal accelerations, gi and g- ln t, decrease. In 
Figure 13, we show snapshots of the surface density of the 
clump and its debris when UMi is on a non-circular orbit 
with apogalacticon and perigalacticon at 89 and 40 kpc, re- 
spectively. For our MOND value, Voo = 170 km s _1 , the 
orbital period of UMi about the Galactic centre is 1.7 Gyr. 
In this simulation, UMi was also initially located at a galac- 
tocentric distance of 76 kpc and has negative velocity, i.e., it 
is approaching to the Galactic centre. The external acceler- 
ation increases until UMi reaches pericentre after 0.45 Gyr. 
An obvious consequence of the adiabatic expansion of UMi 
is that the orbit of the clump becomes more wide with time 
up to 0.45 Gyr. This can be seen clearly in the second panel 
of Fig. 13; the orbital radius is ~ 500 pc at i = 0.5 Gyr. 
By comparing Figs 9 and 13, we find that, at t — 1 Gyr, 



the clump is able to retain more mass in the simulation in 
which UMi is placed on a non-circular orbit about the Milky 
Way but, still, the clump dissipates in < 1.5 Gyr. We can 
safely establish that the varying external acceleration felt 
by the dwarf because of its orbital motion does not offer 
a promising solution to the survival problem of the clump. 
On the contrary, Galactic tides, not included in our simula- 
tions, may help destroy the substructures if UMi's orbit is 
very elongated (say, perigalacticon at < 25 kpc). 

Our simulations have assumed that the mass distribu- 
tion of UMi is spherical whereas the isophlets of UMi have 
cllipticities of 0.54. Since in MOND the gravity comes solely 
from the stellar component, the potential should be natter, 
not spherical. A perfectly spherical potential has orbits that 
are confined to lie on planes. By contrast, orbits in axisym- 
metric potentials generally show precession of their orbital 
planes (except for exactly planar or exactly polar orbits), 
leading to an additional source of orbital mixing. 

We have explored whether a flattened mass distribution, 
such as that observed in UMi, might change the slope of the 
rotation curve in the equatorial plane. In order to estimate 
this effect, we have modeled the stellar bulk of UMi as a 
flattened King model: 



p(R,z) — pK(m) with 



2 n 2 , z 

m -R + — , 



(10) 



where px is the density profile of the spherical King model, 
and derived numerically the circular velocity (see, e.g., §2.3 
in Binney & Tremaine 1987). At radii r < 0.6 kpc, the 
UMi rotation curve in the equatorial plane is again a scaled- 
version of the curve derived in the spherical case. If we adopt 
the same central volume mass density as in the spherical 
case, and a flattening of e = 0.54, the amplitude of the New- 
tonian rotation curve decreases by a factor of 0.85. There- 
fore, if the clump is on a radial orbit in the equatorial plane 
of UMi, the inclusion of the ellipticity is similar to the use of 
a spherical model with Y* = 4.2. Given the uncertainties in 
Y*, this is a secondary effect and a more accurate treatment 
of the flattening is not yet justified by the observations. 

We conclude that the interpretation that the stellar 
clump is a dynamical fossil is difficult to sustain in MOND 
because its lifetime is less than ~ 1.5 Gyr even if clump's 
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Figure 9. Stellar surface density in Mq pc of the debris for a clump on a circular orbit and initial lcr radius of 35 pc and one- 
dimensional velocity dispersion of 1 km s — 1 . Self-gravity of the clump was included. 
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Figure 10. Internal acceleration boost for all the simulated stars 
as a function of the distance to the centre of the clump, at t = 0.5 
Gyr, for the simulation shown in Fig. 9. 
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Figure 11. Distribution of the simulated stars at t = 1.5 Gyr for 
the model shown in Fig. 9. 



self-gravity is taken into account. In the next section, we will 
discuss other scenarios and potential problems. 



6 ALTERNATIVE SCENARIOS AND 
CAVEATS 

In section 5, we consider the dissolution of an unbound 
clump, treating the background as a rigid potential. How- 
ever, dynamical friction may induce a strong orbital decay 
of clusters in dSph galaxies. If the dark matter halo of UMi 
follows a NFW profile, clusters of mass > 5 x 1O 4 M0 sink 
to the centre within one Hubble time if they are initially 
placed within 1 kpc from the dwarf galaxy centre (e.g., Go- 
erdt et al. 2006; Sanchez-Salcedo et al. 2006; Penarrubia et 
al. 2009). In the case of MOND, the orbital decay proceeds 
even faster (Sanchez-Salcedo et al. 2006; Nipoti et al. 2008). 
A fully satisfactory scenario should explain simultaneously 
the survival of the clump against phase-mixing and against 
orbital decay. 

In the dark matter scenario, both problems are solved 
simultaneously if the dark matter halo has a core (Goerdt 
et al. 2006; Sanchez-Salcedo et al. 2006). In MOND, the 
clump can remain in orbit as long as its starting distance is 
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Figure 12. Tidal debris of a clump with apogalacticon at 750 pc 
and on near-radial orbit, at t = 1 Gyr. The clump has an initial 
la radius of 35 pc and one-dimensional velocity dispersion of 1 
km s . 
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Figure 13. Stellar surface density in Mq pc of the debris assuming that the orbit of UMi around the Milky Way has the pericentre 
at 40 kpc and apocentre at 89 kpc. The galactocentric distances of UMi are R gc = (76,41,81,86) kpc, at t = (0,0.5,1.0,1.5) Gyr, 
respectively. In this simulation, the initial clump has a la radius of 35 pc and starts on a circular orbit of 390 pc radius, which is 
represented by a circle in the second panel. Self-gravity of the clump was included. 



large enough so that the dynamical friction force is highly 
suppressed due to the low stellar density of the background 
(Angus & Diaferio 2009). As already discussed, the clump 
cannot survive against tidal diffusion even if its apocentre is 
placed at the tidal radius (~ 2 kpc), because the dissolution 
timescale depends on the gradients of the orbital frequencies 
(e.g., dQ,/dR), which decay slowly with R. 

It is important to consider other alternative scenarios 
for the presence of regions with cold density excess in dSph 
galaxies. It might be that the UMi's clump is not a dy- 
namical fossil but the disrupting remnant of a more massive 
globular cluster in the act of dispersing. Although this possi- 
bility requires a very fine-tuning of the initial parameters of 
the globular cluster and circularity of its orbit, it might be 
viable if the dark matter halo is cuspy because the density 
of dark matter towards the galactic centre might be suffi- 
cient to create a strong tidal field able to disrupt a globular 
cluster of mass < 5 x 1O 4 M0 if it plunges close enough into 
the centre of the dwarf (Penarrubia et al 2009). 

In section 5.2, we saw that a compact clump with ini- 
tial la radius of 12 pc cannot be disrupted by galactic tides. 
Hence, a typical globular cluster of mass ~ 0.8 x 10 5 Mq is 
unaffected by the galactic tides in MOND, even if the orbit 
brings the cluster to the centre of the dwarf galaxy. This 
is expected because the mean density of a typical globular 
cluster is higher than the mean density of a dSph galaxy. 
Therefore, the clump could be a recently disrupted cluster 
which survived for a long time (> 10 Gyr) due to being 
initially bound and more masive, if the parameters of the 
clump progenitor were fine-tuned to be initially much more 
teneous and extended than a globular cluster. How massive 
should the progenitor be in this scenario? The cluster will be 
subject to a prolonged and slow mass-loss regime before it 
suffers from a sharp tidal disruption. Penarrubia et al. (2002, 
2009) and Zhao (2004) find that cored stellar systems lose 
90 — 95 percent of its initial mass during the slow mass- 
loss regime. This implies that the progenitor should contain 
> 13 percent of the total mass of UMi. Clearly, this model 
presents severe inconsistencies because such a massive pro- 
genitor is expected to remain unaltered by tidal stirring. In 



addition, it is difficult to explain why this cluster did not 
sink to the galaxy centre. 

Another related possibility is that the stellar overden- 
sity is not a distinct identity but a portion of a station- 
ary large-scale structure. For instance, density cusps may 
emerge in the projection of eccentric discs (Tremaine 1995) 
or in the projection of the tidal debris of a disrupted stellar 
cluster in triaxial NFW systems (Penarrubia et al. 2009). 
Since simulated cold substructures always show a symmet- 
ric pattern with respect to the dwarf centre, which has not 
been observed so far in the case of UMi's clump, Penarrubia 
et al. (2009) suggest that either a significant fraction of tidal 
debris remains undetected or that the origin of the clump 
is not related to the debris of a stellar cluster. In the par- 
ticular case of UMi's clump, it is unclear if, starting with a 
typical globular cluster, this model can explain simultane- 
ously its high surface brightness and its low internal velocity 
dispersion (Penarrubia et al. 2009). Moreover, the position 
of the resulting stellar substructures in the simulations of 
Penarrubia et al. (2009) are localized at a projection dis- 
tance from the galactic centre of 0.6 — 0.8 times the apoc- 
entre of progenitor's orbit. Therefore, if we assume that the 
orbit of the progenitor in UMi is close to the plane of the 
sky, this constraint places the apocentre of the progenitor 
within two core radius, i.e. ~ 800 pc, which might be insuf- 
ficient to prevent orbital decay of the putative progenitor by 
dynamical friction. 

Wilkinson et al. (2004) suggested that the cold clump 
in UMi is a projection of a cold extratidal population onto 
the face of the dSph. Using numerical simulations, Read et 
al. (2006) explored this scenario in the standard Newtonian 
dark matter scenario and concluded that this hypothesis is 
falsified for two main reasons. First, there is always a source 
of hot tidal stars which mask any cold populations. Secondly, 
the population of stars beyond the tidal radius always ap- 
pears hot in projection because there are more stars on cir- 
cular orbits than radial orbits at the tidal radius. Since we 
cannot see a reason to believe that these arguments cannot 
be applied in MOND, we conclude that if such cold clump 
is real, it could not have formed as a result of tidal effects. 
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7 CONCLUSIONS 

Observations of the rotation curves of spiral galaxies 
strongly support a one-to-one relation between gravity at 
any radius and the enclosed baryonic mass. Due to this em- 
pirical relation, modified gravity theories like MOND are 
able to account rather successfully for the amplitude and the 
shape of the rotation curves. There are some indirect phe- 
nomena that suggest that this empirical relation may break 
down at scales of dSph galaxies (e.g., Gilmore et al. 2007, 
and references therein), challenging the interpretation of a 
modification of gravity as a substitute for dark matter. 

Localized regions with enhanced stellar density and, 
where data permit, extremely cold kinematics have been de- 
tected in some dSph galaxies (e.g., Olszewski &: Aaronson 
1985; K03; Coleman et al. 2004; Walker et al. 2006). In the 
framework of CDM, K03 show that adopting a cored halo 
profile can preserve the UMi's clump incorrupted for a Hub- 
ble time. This cored halo can fit the observed stellar velocity 
dispersion and the persistence of the clump. However, al- 
though we are still unable to make robust predictions about 
how the dark matter distribution changes in the process of 
galaxy formation when the physics of baryons are included, 
the formation of a large core in a dark matter dominated 
galaxy, provides a hard challenge for ACDM. Mashchenko 
et al. (2008) showed that energy feedback in dwarf galaxies 
drives bulk gas motions and gravitational potential fluctu- 
ations large enough to turn the cusp into a flat core (see 
also Governato et al. 2009). They find that for a Fornax 
sized galaxy, the dark matter core has an average density 
of O.2M0 pc~ 3 at redshift z = 5.2, and the density decays 
a factor y/2 in ~ 300 pc. Still, such a core would not be 
enough to avoid orbital spreading of the clump. 

In the particular case of UMi dSph galaxy, the density 
excess around the secondary peak cannot be the remnant 
of a merger with UMi of a smaller, gas-rich system because 
the stars have the same properties in terms of color and 
magnitude as the body of the UMi population (Kleyna et 
al. 1998). One might wonder whether the density peak could 
instead be a projection effect and that what we are seeing 
is a cold, low-density tidal tail. Numerical experiments by 
Read et al. (2006) have shown that this scenario is very 
unlikely. Another possibility is that the clump is a portion 
of the stationary debris of a disrupted globular cluster. This 
model is not implausible but cannot be proven yet. 

A more drastic alternative scenario consists in inter- 
preting the survival of substructures as an internal incon- 
sistency of ACDM. We have explored if MOND can explain 
the longevity of unbound clumps in dSph galaxies. Follow- 
ing K03, we have assumed that this overdensity structure 
is a disrupted stellar cluster and simulate its evolution in 
the gravitational potential derived in MOND. Whichever 
the form of the orbit, the clump is tidally disrupted within 
1.5 Gyr, even if clump's self-gravity is included. One can 
decrease UMi mass by adopting a smaller Y*, slowing the 
process of orbital mixing because the crossing time for the 
clump becomes larger. However, even assuming a mass-to- 
light ratio of O.SMq/Lq, the clump is disrupted in < 2.5 
Gyr. Our conclusion that tidal forces should have disrupted 
the clump appears robust for the adopted value of the stellar 
mass-to-light ratio of UMi. 

The external field acceleration felt by UMi from the 



Milky Way depends on the galactocentric distance of UMi, 
which is time dependent if the orbit of UMi around the 
Galactic halo is elongated. However, the temporal variation 
of the EFE can hardly boost the longevity of the clump. 

In the absence of any alternative model to explain the 
origin of the cold clump, we conclude that it is challenging 
for both ACDM and MOND to explain the nature and dy- 
namics of the clump in UMi. Another alternative of gravity 
suggested by Moffat (2005) is the so-called Modified Grav- 
ity (MOG). It is a fully covariant theory which predicts a 
Yukawa-like modification of Newton's law. For a system with 
a baryonic mass of a few times 1O 6 M0 like UMi dSph galaxy, 
MOG predicts little or no observable deviation from Newto- 
nian gravity at galactic distances of ~ 200 pc (e.g., Moffat &: 
Toth 2008). Therefore, explaining the internal dynamics of 
dSph galaxies is problematic without advocating dark mat- 
ter, weakening the appeal of the MOG model. Other authors 
achieve to find modified theories of gravity which seem to 
reproduce the rotation curves of galaxies (e.g., Capozziello 
et al. 2007). Nevertheless, an analysis of the dynamics of 
dSph galaxies in these theories is still missing. 
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APPENDIX A: MODIFIED EPICYCLES OF 
STARS IN AN UNBOUND SYSTEM IN MOND 

Fellhauer k, Heggie (2005) studied the evolution of an ide- 
alized unbound system in a tidal field. In this Appendix, we 
extend the theory to MOND. The reader is referred to Fell- 
hauer & Heggie (2005) for a discussion about the limitations 
of the theoretical model. 

We consider a small and unbound system on a circu- 
lar orbit with radius Rq in the axisymmetric potential of 
its host galaxy. Therefore, the clump is embedded in the 
external field created by the galaxy $ e . We use rotating, 
cluster-centred coordinates £, C, and z. The £-axis points to 
the anticentre of the galaxy, and the £-axis points in the 
direction of orbital motion of the clump. The presence of 
the clump will change the gravitational potential by some 
increment $;, i.e. $ = $ e + $i. If $; can be treated as a 
perturbation, the gravitational field equation that governs 
the kinematics is 



v 2 + £o|^)<i>, = \- ril 'c;n 



(Al) 



(Milgrom 1986). Here p is the density of the clump, /Uo = 
/x(|V$ e |/flo) and Lq is the logarithmic derivative of po (in 
the unperturbed system). 

We adopt an equilibrium model in which stars are dis- 
tributed uniformly within a triaxial ellipsoid and whose 
epicycles are centered at the centre of the ellipsoid. The 
equations of stellar motion in the epicyclic approximation 
(e.g., Chandrasekhar 1942) are given by 



£ - 20C - 4fiA£ = - 

Z + ** = -ilr> 






Z + V z 



dz 



(A2) 
(A3) 
(A4) 



where £, £ and z are the deviations of the stellar orbit from 
the circular guiding centre. Here v is the vertical frequency, 
v 2 = <9 2 $ e /9z 2 , and A is the Oort constant 



A 



i ( dn 

2\ R dR 



(A5) 



R 



Following Fellhauer & Heggie (2005), it is possible to 
construct a distribution of epicyclic amplitudes so that the 
space density is uniform within a triaxial ellipsoid with semi- 
major axes a, aa and a/Arj, where a > is a free parameter 
that specifies the size of the clump, Ao = VI + Lq, and a will 
be determined from the shape of closed epicycles. The field 
equation (Al) for $i can be transformed into the standard 
Poisson equation by making the substitution £' = f/Arj. The 
resultant equation can be written as 



V 2 ^ = 4^ - 1 Gp(AoC',C^), 



(A6) 



where $• = $■(£', C,z) and V 2 = d 2 /dt, 2 +d 2 /d( 2 +d 2 /dz 2 . 
Therefore, $; is the Newtonian potential created by an el- 
lipsoid of constant density p/po, bounded by the surface 
a 2 = Ao(£' 2 +,z 2 ) + ( 2 /a 2 , which is a prolate ellipsoid with ec- 
centricity e' = \/l — (Aoa)~ 2 . Taking advantage of the the- 
ory of homeoids, $' = 7rG^ ( 7Mip(£' 2 + z 2 ) + nGfi^ 1 pA' 3 ( 2 , 
where A[ = Ai(e') and A' 3 = A 3 (e') are given in Table 2-1 
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of Binney & Tremaine (1987). Therefore, Equations (A2)- is harmonic, we have Q(R) — Oo and n — 2f2, then p cr = 

(A4) can be written in terms of B[ — 2irGpA' 1 /(A 2 ipo) and and k — 0, implying that the dissolution time is infinite. This 

B' A — 2-KGpA^/jj.o as case was already discussed in in terms of the tidal radius in 

i - 2fi< + (- 2 - 4fi 2 )C = -B[£, (A7) SeCti ° n 2 - 

C + 2ne = -^c, (A8) 

z + v 2 z = -A%B[z, (A9) 

where n is the epicyclic frequency if the clump itself is ne- 
glected. As shown by Fellhauer & Heggie (2005), the epicy- 
cle frequency and the vertical frequency are altered by the 
self-gravity of clump: 



11 

K ± = 


i 
2 


and 
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v = 


2 

V 



k +B[+B' 3 ± J(k 2 + B[- B' 3 f + 16Q 2 B! A ,(A10) 



A 2 B[. (All) 

If we choose the lower sign for k' in Eq. (A10), there exist 
exponentially growing solutions (i.e. k' 2 < 0) provided that 

*§<"*■ ^ 

In terms of density, one of the normal frequencies is imagi- 
nary when p is smaller than a certain critical density p cr , 

_ Ag/ip dQ, 2 
P<P -=2^GA^ R 1[R ■ (A13) 

If this condition is fulfilled, the complete solution is a lin- 
ear combination of oscillatory solutions with frequency k! + 
(those obtained by taking the upper sign for k) and expo- 
nentially growing solutions: 

£(t) = £o cos(k+£ + Vo) + A exp(fct), (A14) 

C(t) = -a£o sin(K+t + ip ) - XRcxp(kt), (A15) 

z(t) = zq cos(vt + tpi) , (A16) 



where £o, zo, tpo, tpi and A are free parameters, k > is the 
imaginary part of k'_ , 

2Qk'+ , , s 

(A17) 



and 

R =j?rk- (A18) 

The equations for the epicyclic motions in the Newto- 
nian case are naturally recovered from the above equations 
just by taking po = 1 and Lo = (so that Ao = 1). Due 
to the enhanced self-gravity of the clump by a factor of p^ 1 
in MOND, p CI is smaller than in the equivalent Newtonian 
galaxy, i.e. the galaxy with additional dark matter that has 
the same Q(R). In the particular case that p is much smaller 
than the critical density, we have: 

***(££ -l) B £. (A19) 

We see that the dissolution time for very small initial density 
is oc \/B' 3 oc p,Q 5 p~° in the epicyclic approximation, which 
is valid for small clumps. Since po < 1, the dissolution time is 
shorter in MOND than in its equivalent Newtonian galaxy. If 
the gravitational potential in which the clump is embedded 



